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We study the 5-parameter ("bag factor") for B — B mixing within a recently 
(T) ' developed heavy-light chiral quark model. Non-factorizable contributions in terms 

o ' 

of gluon condensates and chiral corrections are calculated. In addition, we also 

CN 

f-^ , consider 1/mq corrections within heavy quark effective field theory. Perturbative 

fj^ I QCD effects below fx = known from other work are also included. Considering 

OO 



two sets of input parameters, we find that the renormalization invariant i?-parameter 



^ . is B = 1.51 ± 0.09 for B d and B = 1.40 ± 0.16 for B s . 

> 
OO 



O I- INTRODUCTION 

CO 

o 

Studies of the neutral iT-meson system have played a major role in modern particle 
physics [l] . Because of weak interactions, a neutral K meson may be transferred to a neutral 



K meson. This process, known as K — K mixing, determines both the mass-difference 
between the physical neutral states Ki and K$ and the dominating CP-violating effect in 
neutral ii'-meson decays to pions (the e-effect). The neutral S-meson system has rather 
similar properties as the neutral if-system. The difference when going to B — B mixing is 
the importance of other KM quark mixing factors and other mass scales, in particular the 
.B-mesons are about ten times heavier than the i^-mesons. 

In general, non-leptonic processes may be described by an effective Lagrangian which is 
a linear combination of quark operators. The (Wilson) coefficients of the operators can be 
calculated in perturbation theory combined with the renormalization group equations j^. 
At quark level, the leading order diagrams for B — B mixing are given by the so called box 
diagram. This diagram has double W- exchange between two quark lines, and generates 
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an effective Lagrangian(Hamiltonian) for the quark transition bd — > db. This Lagrangian 
has (for all practical purposes) only one operator times a Wilson coefficient containing the 
effects of the virtual (u, c, t) quarks running in the loop. This Wilson coefficient has also 
been corrected for perturbative QCD effects within the renormalization group equations. 
Such calculations has been performed to next to leading order. For B s — B s mixing one 
considers the corresponding bs — > sb transition. 

The difficult part is to calculate the matrix elements of the quark operators between the 
mesonic states, which is a non-perturbative issue. This has been done by lattice simulations 
or by quark models The hadronic matrix element is, as for K — K mixing, 
parameterized through the so called B- ("bag"-) parameter which is by construction equal 
to one in the naive limit when vacuum states are inserted between the quark currents in the 
B — B mixing operator. 

In a previous paper j^, K — K mixing was calculated within a chiral quark model (xQM) 
combined with chiral perturbation theory. Within the xQM, non-factorizable contributions 
can also be calculated in terms of gluon condensates. The purpose of this paper is to perform 
a similar analysis for B — B mixing. We are using a recently developed heavy-light chiral 
quark model (HLxQM) |7[, where non-factorizable effects can be incorporated by means of 
gluon condensates and chiral loops. 



II. B - B MIXING AND HEAVY QUARK EFFECTIVE THEORY 



At quark level, the standard effective Lagrangian describing B — B mixing is 
2 = - -§ M hV (Y: b V tq f S (x t ) tjb 6(a0 Q(AB = 2) , 



4vr 2 



(1) 



(for which q = d or 
due to short distance 



where Gf is Fermi's coupling constant, the Vs are KM factors 
s for Bd and B s respectively) and Sq is the Inami-Lim function 
electroweak loop effects for the box diagram: 

4x — llx 2 + x 3 3a; 3 Logs 
4(1 -x) 2 2(1 -xf ' 

In our case, x = x t = mf/M^, where m t is the top quark mass. Because of its large mass, 
the top quark gives the dominant contribution. Also the u and c quarks are running in the 
loop, but these contributions are KM suppressed. The quantity Q(AB = 2) is a four quark 



(2) 
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operator : 

Q(AB = 2) = ql 1 a b L ql la b L , (3) 

where qt is the left-handed projection of the q (fr)-quark field. The quantities T] B 
and b(n) are calculated in perturbative quantum chromodynamics (QCD). At the next to 
leading order (NLO) analysis it is found that r] B — 0.55 ± 0.01 |2j. Furthermore, for a 
renormalization point /z in perturbative QCD equal to or below m^, 



&0) = Mi*)] 6/23 



1 + ^T Js 



(4) 



where J5 = 1.63 in the naive dimension regularization scheme (NDR). At fi = nib (= 
4.8 GeV) one has &(m&) ~ 1.56. 

The matrix element of the operator Q(AB = 2) between the meson states is parameter- 
ized by the bag parameter B Bq : 

(B\Q(AB = 2)\B) = lf B M 2 B B Bq (fi) . (5) 

By definition, B Bq = 1 within naive factorization, also named vacuum saturation approach 
(VSA). This means to insert a vacuum state between the two heavy-light currents in the 
operator Q(AB = 2), and use the matrix elements defining the decay constant j B : 

(0\qE^b\B(p))= l -f B p» and (B(jp)\qE^b\0) = -- f B p^ . (6) 

One may combine naive factorization with the large N c expansion, where A^ c is the number 
of colours. Then one finds B B = 3(1 + l/iV c )/4, giving B Bq = 3/4 in the (naive) large 
N c limit. We will see later that there are important non-factorizable contributions of order 
1/N C . In general, the matrix elements of the operator Q(AB = 2) are dependent on the 
renormalization scale /1, and thereby B Bq depends on [i. As for K — K mixing, one defines 
a renormalization scale independent quantity 

B Bq = b{iS)B B U) . (7) 



Within lattice gauge theory, values for B Bq between 1.3 and 1.5 are obtained 

The mass difference between the weak eigenstates (B H and B L ) are related to the bag 
parameter in the following way for B q = Bd,B s : 

Am o = ^™B q fl q B BqVB M 2 w So (ml/Mly) \V: q V tb \ 2 . (8) 
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In order to extract the KM matrix elements it is crucial to have a precise knowledge of the 
bag parameter B Bq , and the weak decay constant fs ■ 

The b -quark is heavy compared to the typical hadronic scale of order 1 GeV, where 
confinement and chiral symmetry breaking effects are essential. Perturbative effects below 
the 6-quark scale may then be calculated down to f GeV by means of heavy quark effective 
theory (HQEFT. See Q| for a review). Thus HQEFT also allows us to evolve the matrix 
element (J3J) from \x = rrib down to 1 GeV. 

HQEFT is a systematic expansion in l/m&. The heavy quark field b(x) is replaced by a 
"reduced" field, Qv (x) or Qi \x), which is related to the full field the in following way: 

Q&(x) =P ± e™- x b(x), (9) 

where P± are projecting operators P± = (1 ± 7 • v)/2. The reduced field Qv can only 
annihilate heavy quarks. In order to describe heavy anti-quarks one has to use Qv } . In 
other words, Q^\Q v ■*) annihilates (creates) a heavy quark (anti-quark) with velocity v. 
The Lagrangian for heavy quarks is (Q v = Qv)' 

where _D M is the covariant derivative containing the gluon field (eventually also the photon 
field), and a ■ G = o^ v G a ^ u t a , where a^ u = i[y^, 7"]/2, is the gluonic field tensors, and t a 
are the colour matrices. This chromo-magnetic term has a factor Cm which is one at tree 
level, but slightly modified by perturbative QCD effects below the scale m h . It has been 
calculated to NLO Furthermore, {iD ± ) 2 eS = C D {iDf - C K {iv ■ Df. At tree level, 

Cd = Ck = 1- Here, Cp is not modified by perturbative QCD, while Ck is different from 
one due to perturbative QCD corrections |U|. In our case, vhq = rrib is the heavy quark 
mass. 

Running from [i = mi, down to /i = A x = I GeV, there will appear more operators. Some 
stem from the heavy quark expansion itself and some are generated by pertur 
effects. The AB = 2 operator in equation (J3J) for A x < \i < m h can be written 

Q(AB = 2)= C 1 (jM)Q 1 + C 2 (jM)Q 2 + —(y ( a i (ji)S i (jM) 

+ £ hiMXiQi) ) +0(l/m 2 b ) . (11) 



£ H qeft = ±Q V tv-DQ v + ^-Q v (-C M |^ ■ G + (iD ± ) 2 eS j Q v + 0(m Q 2 ) , (10) 



jative QCD 



14. 



15|,lig : 
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The operator Qi is Q(AB = 2) for b replaced by while Q 2 is generated within pertur- 
bative QCD for fi < nib- The operators Si and Xi are taking care of 1/mj, corrections. The 
quantities Ci, C 2 , Oj, hi are Wilson coefficients. (C\ = 1 + 0(a s ) and C 2 = + 0(a s )). The 
explicit expressions for the operators are 

Q 1 =2qEYQl +) Wl,Q { - ] , (12) 
Q 2 =2gZ^gW qZv,Q ( f ] , (13) 

-2ie Xfiup v x qEiD^Q[ +) 1E1 p QV 

+ 2qZQ[ +) qZiQQV + 2gZ7M^l +) W^Q^ 

-i2e Xfivp v x qZYQl +) qEiD^Y , (14) 
X 2 = 8[iv d(wQl +) )] qZQV +2[iv d(Wl»Ql +) )] WYQ[' ] , (15) 
X 3 = 4 [iv ■ d{q- L l, Ql +) )] qll" QV • (16) 

The operators Si are nonlocal and is a combination of the leading order operators Qi,2 and 
a term of order I/uiq from the effective Lagrangian (|10|): 

5l • ^d^QiCO),©^)}, 

d 4 yT{Q 2 (0),O K (y)} , 

d 4 yT{Q 1 (0),O M (y)} , 

d 4 yT{Q 2 (0),O M (y)} , (17) 



m 6 

^2 _ , / .,! 



m 6 



where 
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^ a • GQ[ +) + Q< -) a • G?QH) , (18) 

are the kinetic and magnetic operators of eq. (JTDJ). There are no mixing between the local 
operators and the non-local operators, since the local operators do not need the non-local 
ones as counter-terms. The Wilson coefficients aj will then be the product of C\ 2 and Cm k- 



The Wilson eocfneionts ( \ and ( ' 2 have been calculated to NLO [1J, Q| and for fj, = A x , 
Ci(A x ) = 1.22 and C 2 (A X ) = —0.15. The coefficients ^1,2,3 have been calculated to leading 
order (LO) in [l^, and the result at \i — A x is h\ = 0.52, h 2 = —0.16 and h 3 = —0.15. 
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III. THE HEAVY-LIGHT CHIRAL QUARK MODEL 

In order to calculate the matrix elements we will use the heavy-light chiral quark model 
(HL^QM) recently developed in Q|. This is a type of quark loop model jl^, 13, Q, Q 
where the quarks couples directly to the mesons at the scale of chiral symmetry breaking 
A x , which we put equal to 1 GeV. What makes our model [7j distinct from other similar 
models is that it incorporates soft gluon effects in terms of the gluon condensate with lowest 



dimension 



21 



22. 2>lQ- 

Lagrangian describing this interaction can 
be obtained as a mean-field approximation of the (extended) Nambu-Jona-Lasinio model 

(njl) 0,0- 

In this section we will give a short presentation of the HLxQM. In the next section we 
will use the model [7] to calculate non-factorizable soft gluon effects in B — B mixing. 



The Lagrangian for the HL^QM is 



C 



HL X QM 



C 



HQEFT 



£*QM + Ant 



(19) 



The first term is given in equation ([10)1 . The light quark sector is described by the chi- 
ral quark model (xQM), having a standard QCD term and a term describing interactions 
between quarks and (Goldstone) mesons: 



C 



xQM 



X [f(iDn + V M + -m} X -xM qX 



(20) 



where Xl,r are the flavour rotated quark fields given by: 



Xl = tfqL 



XR = &R 



(21) 



where q T = (u, d, s) are the light quark fields. The left- and right-handed projections 
and qn are transforming after SU(3)l and SU(3)r respectively. The quantity £ is a 3 by 3 
matrix containing the (would be) Goldstone octet (tt, K, r/) : 



£ = e m// where IT = — a (x) = 

2 v 2 



K- 



7T 



— EL _i_ HI. 

V2 ^ s/E 



K+ 
K° 



--4 



(22) 



where / is the bare pion decay constant. In 1)20)1 . m is the (SU(3) - invariant) constituent 
quark mass for light quarks, and M q contains the current quark mass matrix M. q and the 
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field f : 



M g = M v q + Mf l5 , where 

1 



(23) 



Mf = -@M$+ ZM£) and Mf=--@Mtf-ZM£). (24) 



The vector and axial vector fields V M and in (|2()jl are given by: 



(25) 



Furthermore, the covariant derivative -D M in (|2()|) contains the soft gluon field forming the 
gluon condensates. The gluon condensate contributions are calculated by Feynman diagram 
techniques as in 22l |23J] . They may also be calculated by means of heat kernel techniques 



as m 



21 



25 



The interaction between heavy meson fields and heavy quarks are described by the fol- 
lowing Lagrangian : 

1 



C 



Int 



-G 



H 



H^Q^+Qt ) H^ Xa 



+ 



2G; 



-Tr \H*Hfi 



(26) 



where Gh and G3 are coupling constants and is the heavy meson field containing a 
spin zero and spin one boson: 



(27) 



The fields p( + \p(~^) annihilates (creates) a heavy meson containing a heavy quark (anti 
quark) with velocity v. 

Integrating out the quarks by using (fTUj) . (j2*U|l and (pjj) . the effective Lagrangian up to 
O^q 1 ) can be written as 



ixks by " 



C = TTr 



Hp* (™-D 6o -A Q )# 6 (±) 



QaTt 



Tt{±)tt(±) A,! 



(28) 



where iD^ a = i5b a D^ — V^ a . The term proportional to the quark- meson mass difference 
Aq = Mh — ttlq in ()28|) is irrelevant for us due to the reparametrization invariance jlo| . 
Also, it does not enter our loop integrals because our heavy meson fields are attached to 
our quark loops at zero external momentum. (The external momentum includes the piece 
v^Aq). As shown in B, the term ~ I/G3 in (J26j) is related to Aq, and this term is also 
irrelevant within the present paper. 
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To obtain (|2g|) from the HL%QM one encounters divergent loop integrals, which might 
be quadratic-, linear- and logarithmic divergent. For the kinetic term in ()28|) we obtain the 
identification: 

- iG%Ne ^ /2 + 2ml2 - i ^. { ^))^, (29) 

where J3/2 and I 2 are the linear and logarithmic divergent integrals respectively, and {^-G 2 } 
is the gluon condensate. To obtain the axial vector term proportional to g A , we obtain a 
similar condition, and combining it with (ffijj) . we obtain for the axial vector term 

9A = l + ^G 2 H N c (l 3/2 -^\ , (30) 

such that the (formally) linear divergent integral 1 3/2 is related to the strong axial coupling gj± 
(or strictly speaking, its deviation from one). Analogously, within the pure light quark sector 
(the xQM), it is well known that the quadratic and logarithmic divergent i nteg rals are related 
to the quark condensate and the bare decay constant /, respectively 12, 21, 2^,0, 

(qq) = -4imN c h - — ^— ( — G 2 ) , (31) 

12m 7r 

P =-^NJ 2 + ^- 2 (^G*) . (32) 

The divergent integrals Ji, J 2 and / 3 / 2 are listed in appendix |XJ The effective coupling 
Gh describing the interaction between the quarks and heavy mesons can be expressed in 
terms of m, /, g^, and the mass splitting between the 1~ state and 0~ state. Using (|2~H|) . 
(J3U|) . (|32j) one finds a relation between this mass-splitting and the gluon condensate via the 
chromomagnetic interaction in (jlU[) j?| : 



W 2 /4 ^ 2 _2m„ m _ (vr + 2) 



where 



( V G2) = ^T' G - = 7^' ^^(A x)) (33) 

(1 + 3g A ) + 3 
P = 4n ! ' &{H) = U Q {M H * - M H ). (34) 

In the limit where only the leading logarithmic integral I2 is kept we obtain : 

9a - 1 , P - 1 , G H - Gg } = ^ . (35) 



Note that g A — 1 is the non-relativistic value [27]. We observe that the mass-splitting 
between H and if* sets the scale of the gluon condensate. This means that, while in |23j 
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the gluon condensate was fitted to the K — > (27i) I=2 amplitude, it is here determined in the 
strong sector alone (with a slightly lower value than in 



n 

The 1/rriQ corrections to the strong Lagrangian have been calculated in pj. They may 
formally be put into spin dependent renormalization factors. This means that (|28jl is still 
valid with the replacement H T = H(Zh)~z, where Zjj and the renormalized (effective) 
coupling «7_4 are defined as: 

Z H = 1 + , 36 

m Q 

9 A = 9a(i — (ei - 2d^e 2 )j (#i - d A g 2 ) , (37) 

V m Q J m Q 



where 



(38) 



and 



3 for 0" 1 for H*H coupling 

i.m \ d A = < 

•1 for T —1 for H*H* coupling 



" 1 4m J 16vr 16 V m J ' 

+ 12^ (C - + 8 -^<V G2 >) • < 39) 

9l = m - d, | m + 1 + N ^ + *) _ ql ( <»> + 3f) 

'' 1 12m 6 48tt 16 m J ' 

' -(C i ,-2 7 r)(^G 2 )') , (40) 



64m 2 7r 

<7 2 = 7^4^ • ' (41) 

y2 (7r + 2)6m V ; 

IV. BOSONIZING Q(AB = 2) 

In this section we will discard I/tuq terms. We are then left with the operators Qi,2 
defined in equation (|12j) and (J13)) . In order to find the matrix element of Q12 , one uses the 
following relation between the generators of SU(3) C I, n are colour indices running from 
1 to 3): 



SijSin = —Sinhj + 2 t a in tfj , (42) 



1 

N c 

where a is an index running over the eight gluon charges. This means that by means of a 
Fierz transformation, the operator Q\ in (|12J) may also be written in the following way : 

Qi = ^WYQ^Wi.Q^ + mt a i> x Q i + ) qZt a ^Qi- ) , (43) 
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FIG. 1: Diagrams for bosonization of the left handed quark current 
and similarly for Q 2 . 

The first (naive) step to calculate the matrix element of a four quark operator like Qi 
is by inserting vacuum states between the two currents. This vacuum insertion approach 
(VSA) corresponds to bosonizing the two currents in Qi and multiply them, as mentioned 
below eq. (JSj). For one current, visualized in figure ^ one obtains 



tiply the 
AH : 



U-fQ^ — °^Tr [d f l a LHt } \ , (44) 

Using the relations - for the divergent integrals, and also eq. (JHB1), we obtain 0: 

Gh ( (qq) -,2/1 1 \ , ( 7^ - 2 )/ a ^2^ 



clh = ^ " 2f (1 - -) + ^r^<-^> • (45) 

This bosonization has to be compared with the matrix elements defining the meson decay 
constant f B given in eq. ©• In those relations, b is the full quark field. Within HQEFT 
this matrix element will, below the renormalization scale fi = rriQ (= mj,), be modified in 
the following way: 

(0\qZ^Q^\B(p)) = % -f B M B v» and (B(p)\qZT» Q^\Q) = ~f B M B v" , (46) 
where Q 

= CM^ + C v (fi)v^ . (47) 
The coefficients C 7i „(/i) are determined by QCD renormalization for /z < ran . They have 



been calculated to NLO and the result is the same in MS and MS scheme |28l |. In FILxQM 
the decay constant f B can be calculated and the result is |7[ : 

aH c^ + cm g» ' 1 J 

The second matrix element in (J4H)) is genuinely non-factorizable, and we have to go beyond 
the VSA. However, in the approximation where only the lowest gluon condensate is taken 
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into account, the last term in (f4*H|) can be written in a quasi-factorizable way by bosonizating 
the heavy-light coloured current with an extra colour matrix t a inserted and with an extra 
gluon emitted as shown in figure [21 Calculation of this diagram is straightforward when 
using the light quark propagator with just one soft gluon emitted : 



S G (k) = jG b a / [<r a/3 ( 7 • k + m) + (7 • k + m)a a P] (k 2 - m 2 )~ 2 . (49) 



9s n b + b \„aBi 

The part of the diagram to the left in figure 121 then gives the bosonized coloured current: 

, (50) 



{Wt a l a Ql ±] ) 1G — ► -^G^Tr ^LH^UiI 2 {a^,j.v}+ ±-oA 



where I2 is to be identified with f 2 by the use of equation (JH2*)) . The result for the right part 
of the diagram with B replaced by B is obtained by just changing the sign of v and letting 
P 5 — > Pg (remembering that Pg ^ creates a meson with a heavy anti quark). Multiplying 
the coloured currents, we obtain for the non-factorizable part of Q\ and Q2 to first order in 
the gluon condensate : 



Ci Wt a YQi +) Wt a 7„ Qi~ ] + C 2 Wt a qlt a v, QH 



where 

^2 f / -f \ 2 o 2 



128 I N c \mJ N? V m 



and S = £ 2 , where £ is given in Eq. (|22|) . Note there is no sum over i, i = 2,3 for q = d,s 
respectively. 

The Lagrangian in equation (|2*UJl contains couplings involving the the current mass term 
and the chiral quark fields. This makes it possible to calculate the counter-terms needed in 
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FIG. 3: Mass insertion in the nonfactorizable part of the current 



oops. The counter- 

r 

3] when calculating 



order to keep the chiral Lagrangian finite after the inclusion of chiral 
term for the factorizable part of the amplitude has been considered in 
fs- In the case of the non-factorizable part of the amplitude, we need to consider similar 
diagrams as those shown in figure |21 with mass insertion like in figure where mass insertion 
is indicated by a cross on the light quark line. The bosonized current with mass insertion is 



32mvr 2 ^ - ' - i ">/ i > < ■ ( 53 ) 

This result can also be obtained by simply differentiating the right hand side of equation 
(jnOJ) with respect to m. 

The bosonized version of the Q(AB = 2) operator can then be split in a pseudo scalar 
and a vector part: 

Q(AB = 2) Bos . = ApP^TbP™ + AyPt^tA^ , where : 

A ? = li 1 + W)( C 1 ~ C *) a H f 1 + 2 — m «) ~ Ci(-G 2 ) ((3 B + u p m q ) , (54) 
2 JM C \ a H J 7r 

A v = hl + ^)C ia 2 H (l + 2^m q ) - (-G 2 ) ({C, - %fi B + C W 
2 JM C \ a H J 7r V 6 

The quantity up is the counter-term obtained from (J55)l . and uj\ is a counter-term for f Bs 

found in Q: 

G 2 H f 4tt/ 2 1 . . 

^ = 6^\ 1 + iW/' (55) 

(1-3&0 (9tt - 1Q)G H a s 2 
^ = —G~ H T92^ ( V G > ' (56) 

For the current quark mass entering (JH3|) we will use 

m d = -mlf 2 /(qq) , and m s = -m 2 K f/{qq) . (57) 
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The term including the vector fields P M are needed in order to calculate chiral corrections 
where B* are included. From equation the equations ©, © and (jSljl the renormalization 
invariant bag parameter can be extracted. Anticipating the results of the two next sections, 
it can be written in the form: 



B 



B 



3 -l 
4 



1 



where 




+ — + ( I 



N C J 32n 2 f 2 



We find from (|54|) the parameter due to genuine non-factorizable effects: 



() G 



7T 



a 



H 



2d 



C\ — C*2 



(58) 



(59) 



(60) 



Note that this parameter is formally of order (N c )° and is positive, which means that this 
non-factorizable contribution reduces the value of B according to (|58j). Thus we are quali- 
tatively in agreement with |5[, where a negative contribution to the bag factor from gluon 
condensate effects is found. 

Using the relation between an and fs in Eq. ()48|) and the expression value for Gh in 
equation we may also write : 



SB 




32tt 2 / 2 /1 M b 



Ci 



C\ — C*2 



(61) 



Numerically, / and fs are of the same order of magnitude, and 5q is therefore suppressed 
like tu/Mb compared to the corresponding quantity 



Nr. 



; 32tt 2 / 4 ' (62) 
for K — K mixing. However, one should note that fs scales as l/yMe within HQEFT, and 
therefore 5q is still formally of order (mi,) . 

The formula (|58|) is a generalization of a similar formula found for K — K mixing 6]. The 
quantities r& and will be calculated in the next sections, while r x is known from previous 
work (2^. More specific, the quantity r^, to be calculated in the next section, has dimension 
(mass) 1 and depend on hadronic parameters calculated within the HLxQM. Similarly, the 
quantity r x contains the chiral corrections to the bosonized versions of to be presented 
in section VI. The quantity contains the chiral corrections proportional to {^-G 2 ) and 
the counter-terms up and iO\. 
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V. 1/mg CORRECTIONS 

The 1/mg corrections have been defined in equation (jl4H17j) . In the HL^QM we only 
need to consider ()14j) and (fT7j) . This is due to the fact that when we are considering terms 
in the effective Lagrangian for B — B mixing the external particles carry no redundant 
momenta |7[. (In other words, the 5-mesom momenta are ps = Mbv). Hence the operators 
in (fTo*)) and ()16|) will give zero contribution. 

The operator in equation (fT4*|) can be written on the form 

3 6 

X 1 = 2^qT j iD a Q<f)qe j QP + 2 £ qTj Q« q 0, iD a QH , (63) 

3=1 i=4 

where r Q , are defined : 





= RY 




0i 


= R 


r 2 


= Rg^ a 




©2 


= Ri» 


r 3 


= -i e Xaup v x 


RY 


©3 


= RY 


r 4 


= R 




©4 


= RY 


r 5 


= R 




©5 


= Rg^ 


r 6 


= -i e Xav P v x 


RY 


©6 


= Rj u , 



(64) 

where D a is the covariant derivative containing the gluon field. Note that the operator Xi 
is Fierz symmetric Q|. We bosonize X\ in the same way as Q\.2- 

Some two-quark operators appearing in (|B^jl are already studied in \\ when calculating 
l/m b corrections to We use those results when bosonizing X 1; and the result can be 
written: 



X\ - x\ 

3 



bos 



£ {2(1 + h^-Tr [^QiH^] l -Tr [^H™ (aJY + 
i=i < c 

+ 4f3.Tr [tfQiHU (-& {Y\ 7 ■ v} + faY") } Tr [^H^ ((3 3 D, ua + 2mfca^v a ) } } 

+ £ W l + w )2 f TT [$ V * H P («37 Q - irr [e f ©^i +) ] (65) 

i=4 ^ c 

+ 4&Tr [^r^H (&£W - 2mftv« a )] Tr (ft K", 7 ■ ^ } + ft<0] ) } , 

where = {cr M „, 7^} (g^ — v a vp). The second and fourth lines are genuinely non- 

factorizable. The a's and /3's are hadronic parameters calculated within the HLxQM, and 
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are given in Appendix [BJ Evaluating the sums and traces in equation (|65p we arrive at : 



bos 



(2) 

where ^ is a combination of the /3j's and can be written 



7T 



-P 



(~>yt p(+) 

a in 



Qp(-)yt p(+) 



7T 



1+ 



4tt // 



(66) 



(67) 



3N C \m. 

The bosonating of the nonlocal operators is rather straight forward in this model. The 
result for the factorizable part of the non local operators can be found in [t| in the calculation 
of/fl : 



E 



i=i 



AiSf 8 * 
m b 



- 1 



1 



N C J m b G H 



//m /~f \ p( — )v*t p(+) 



(68) 



The result for the nonfactorizable part of the operators is : 



E 

i=i 



A cNfact 



— (—G 2 )p K ( (d - ic7 2 )P/-^t P« + CxPHElpW) 



M 



1 a 

where the quantities /3# and s are given in Appendix iBl 
We need fs which has been calculated in [tJ to l/m& : 

K 



p(->yt p(+) _|_ qp(-M p(+) 



(69) 



V m 6 



«6 



32tt 2 / 2 

(ei-6e 2 ) , (S^ + ^aS) (/4 ~ /4) 



where : 



2a#(C 7 — C„) 



GnotH 



"Xd 



"Xs 



m 



-m l K {\ + ^ 2 ) + m 2 K (ln + In ^)(1 + 3g A 



11 



11 
18 
13 

— { -m 2 K (l + g A 2 ) + m 2 K (ln + ^ In ^)(1 + 3^ 



(70) 



m|- 4.4, 



, ^32tt 2 / 2 
H m s 



^ 13 3' 

where P 7 and P^ are sums of Wilson coefficients. The contribution to the bag parameter 
from l/m& corrections can now be extracted (see eq. (jKHjl): 



1 + 



6Ci 



6Pi 



P, 



a# \Ci — C2 C 7 + C t 

2\ J ^1 o(2) , 



Pj; 



(Ci - C 2 )a 



7T 



G 2 ) -±fi> + ^ + CmA 



(i) 

M 



G x 1 ■ 



M 



(71) 



It should be noted that l/rrib corrections increases P, in agreement with [l 
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FIG. 4: Diagrams contributing to the bag parameter 



VI. CHIRAL CORRECTIONS 



We will only consider chiral corrections to Qi^ in equation (|T2"|) and ()13jl. Adding chiral 
corrections to operators proportional to 1/tuq will be considered as higher order. The chiral 
corrections to the bag parameter have been considered in 2^|. Some of the corrections are 
simply corrections to jsfll. Isil. 8^ |. The diagrams shown in figure are those which are 
genuinely non-factorizable, i.e. they are not included in chiral corrections to /# . 

The chiral corrections (r x ) to the bag parameter can then be written : 



(L 



9 
Ci 



m K In 



Ci-C. 



-0A 



Am 2 K 
-m 



K 



r, 



G 



(L 



9 A 



4m|- 



9 K 
A 2 ) In 

2 o 



"V" 



--m 2 Y + -A 2 (2-3F(A/m r , 



-m 



9 



K 



Ci - C 2 /3 



-gL 



&4 



, +2 



-m^ - A 2 ) In 

327r 2 / 2 m s , 



4m 



K 



3/i 2 



-m 



-A 2 (2 - 3F(A/m v ) 



(72) 



where we have ignored the pion mass and used the mass relations m 2 



(73) 

4m 2 K /3. The 
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T 



#(+) J ' M 



2ig A k fJ , Ili 



T 



2ig A k^ U i3 



H{k) 



i • = 



£P(&) 

i, M • = = = = =* J> 



2(vk- A) 



(5,., 



' 2(wfc-A) ' 



FIG. 5: Feynman rules for the strong sector, II is given in equation 1)22(1 



function is denned in equation (|A7|) and : 



1 for B rl 



4 for 5, 



and gL 



for 5 d 

1 for B, 



(74) 



If one ignores the counter-ter m g iven by up, and take the limit A = M* H — Mh — > 0, we 
obtain the same result as in [23. For the bare coupling constant / we will use the value 

rn 

/=86 MeV |32j. The Feynman rules for chiral loops are given in figure El 



VII. NUMERICAL RESULTS 



n 



The model dependent parameters of the HLxQM was fixed in |7J by using various 
constraints. For instant, the constituent light quark mass was determined to be m = 
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FIG. 6: The bag parameter B for 



0.28 



0.30 



220±30 MeV. Using the parameters from []|, 

n = (0.26 ± 0.04) GeV 
-(0.02 ±0.01) GeV 2 
-(0.03 ±0.01) GeV 2 
= 1.53 ±0.05 
f Bd = (170 ±25) MeV 



T- 



Ari 



we obtain (using A = M* H -M H = 0.025 GeV): 

(0.5 ±0.1) 
-(0.10 ±0.04) GeV 2 



SB 



T: 



T. 



G 



t° = (0.12 ±0.06) GeV 2 



Yd 



X 

Bb 



B 3 



1.48 ±0.08 
(180 ± 25) MeV 



f Bd jB B , = (215 ±30) MeV 



f Bs - 



SB d \/B Bd 



1.05 ±0.01 



f Bs \/B Bs = (225 ±30) MeV 



1.08 ±0.02 



(75) 



The decay constants f Bd and f Bs were also given in 0, but are listed also here for com- 
pleteness. (Note, however, that the values are slightly different, because in [3| we did not 
distinguish f n from the bare coupling /.) The values for the bag parameter B are in agree- 
ment with lattice calculations 0, Q. A plot of B as a function of the constituent quark 
mass m is shown in figure |H1 and [7| We observe that the values of B are fairly stable over a 
large variation of light quark constituent mass m. Especially this is the case for B^. From 
m = 180 MeV and m = 300 MeV the bag factors only changes with 10%. We note that 
1/rrib corrections are small. 

The values for the f B 's and especially for the ratio f B J f Bd (and £) in (|73j) are a bit 
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low 0, 0|. There might be at least three reasons for this. First, concerning the absolute 
value for /#'s, they dependent significantly on the value of the quark condensate, as seen 
from equation (jEJ) and In Q we used the "standard" value (gg) = (-240 MeV) 3 , 

without any uncertainty. It could be argued that we should have used an uncertainty 

1 /3 

of 10 MeV, say, for (gg) , although the wide range 190 to 250 MeV used for m will 
to some extent compensate for this. Second, it might be that our expansion within the 
HLxQM overestimates the counter-term U\ which reduces f Bs . However, neglecting this 
counter-term would give the high value fnJfB d — 1-3. Third, our value for the axial pion 
coupling gj± in (f2~£J) might be too low. In we used input from QCD sum rules Q] both 
in the B- and .D-sectors. Alternatively, we may use the experimental value for the effective 
coupling g_A H * Hlx = 0.59 ± 0.09 in the D-sector j^, giving almost the same bare coupling 
g_A = 0.59 ± 0.04. Using this bare coupling also in the 5-sector (instead of g_A = 0.42 ± 0.06 
in 0]), and in addition (qq) 1 ^ 3 = (—240 ± 10) MeV, we obtain an alternative set of values: 

r b = (0.25 ± 0.04) GeV 5 G = (0.5 ± 0.2) 

r Xd = -(0.06 ± 0.01) GeV 2 r Xa = -(0.25 ± 0.04) GeV 2 



T ^ d = -(0.07 ± 0.01) GeV r£ = (0.2 ± 0.2) GeV' 

B Bd = 1.51 ± 0.09 B Bs = 1.37 ±0.14 

f Bd = (190 ± 50) MeV f Ba = (210 ± 70) MeV 



f B \JB Bd = (240 ± 70) MeV f Bs \j B Bs = (260 ± 90) MeV 
f = = 1.08 ±0.07 ^1 = 1.14 ±0.07 (76) 



We observe that the value for f B J f Bd in (fTo^l is close to the standard one. 

To conclude, we have calculated the bag parameter B for the B^ and B s mesons. 
Combining our two alternative sets of values (and consider the range of values) we find 
B Bd = 1.51 ±0.09 and B Bg = 1.40 ±0.16. The value for B Bs is more sensitive to chiral loops 
and counter-terms, and therefore the uncertainty is bigger. 

In principle, B is renormalization invariant (/i independent). This cannot be shown 
within our approach. By construction, perturbative QCD within HQEFT, the HL%QM and 
chiral perturbation theory are matched at the scale A x . However, we have a reasonable good 
matching numerically as in (23J. Varying the renormalization scale /i = A x in the range 0.8 
GeV to 1 GeV, the bag parameters only change with 6%. Moreover, like in p, the formula 



2.00 
1.80 



1.20 
1.00 



1 


1 1 

LO - 


- 1/mQ + x - 








LO + x 




- 


Leading 


Urder (LU ) — 


- 


1 


1 1 


1 





0.18 0.20 0.22 0.24 0.26 0.28 0.30 



to (GeV) 



FIG. 7: The bag parameter B for B. 



(I58J) nicely shows the the various parts building up the total result for B. 



APPENDIX A: LOOP INTEGRALS 



The divergent integrals entering in the bosonization of the HLxQM are defined : 

, d d k 1 

h 
h 



(27r) d k 2 - m? 
d d k 1 
{2Ti) d (v ■ k)(k 2 -m 2 ) 
d d k 1 
(27r) d (k 2 - m 2 ) 2 

The integrals needed in the calculation of chiral corrections to the bag parameter are 

a f d d k 1 — zA (\ , . 9 , . i k \ 

L™f = / - — -— — — = — - - - n (m 2 +2 -2F m/A 

1,1 J (2n) d (k 2 - m 2 ){v ■ k - A) 8vr 2 \e V 1 V 1 1 
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Ag^ + BvW 



A 



(2ix) d (k 2 - m 2 )(v ■ k - A) 

1 f d d k k 2 -{vk) 



d-lj (2n) d (k 2 - m 2 )(v k - A) 



/A \-- + ln(m 2 ) - l)(m 2 - -A) - -F(m/A)(A 2 - m 2 ) - -(m 2 - -A 2 ) XA5) 



16tt 2 I v f v " ' /v " 3 7 3" v "' /v "" y 3 V " 6 

, d d k (vk) 
B = -A+ ' 



(27r) d (k 2 -m 2 )(vk - A) 

Zl^ ((_! + i n ( m 2 } _ 1)(2m 2 _ 8 A) _ 4 F(m/A)(4A 2 _ m 2, 
16n z I e 3 3 

-^(m 2 - ^A 2 ) \ (A6) 



where: 



— y/x 2 — 1 tan 1 (\/x 2 — 1) x > 1 

/•'(■'■)= < (A7) 

Vl - x 2 tanh - x 2 ) x < 1 



In the case of A > m we have ignored an analytic real part in (jA4|) . Equation ()A4|) coincides 



with the one obtained in j^J however equation ()A6)) differs by a factor — 2/3 (m 2 — 2/3 A 
inside the parenthesis of the expressions for A and B. This is presumably due to the factor 
l/(d- 1) = (l-2/3e)/3 in A. 



2\ 



APPENDIX B: SOME DETAILED EXPRESSIONS FOR HADRONIC 

PARAMETERS 



The parameters of equation (J55J) are : 



m G H _ 

« 3 =-3 + ~g~\99/ 

777- 2 

«3 = + ^G H {qq) 

ft 

12 7T 



The /3^'^'s in (pj§j) are given by : 



f 2 



Am 2 N c 
AG 2 B N C 



(Bl) 
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_m_ 2 I 4_8tt //\ 2 1 32tt 2 / ' x 1 

256^1 1 + tt j\£ I m J 1 p ^ iVf V"m 



r J M 



3 (2) 



7T 



12N 2 U 



2 



87T m' + 167T 2 



N c \m ) iV 2 \ m 



7T 



24N r 



1 + 



2tt // 



N c \m 



(B2) 

(B3) 
(B4) 
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